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In fluid transport across nanopores, there is a fundamental dissipation that arises from the connection between 
the pore and the macroscopic reservoirs. This entrance effect can hinder the whole transport in certain 
situations, for short pores and/or highly slipping channels. In this paper, we explore the hydrodynamic 
permeability of hourglass shape nanopores using molecular dynamics (MD) simulations, with the central pore 
size ranging from several nanometers down to a few Angstroms. Surprisingly, we find a very good agreement 
between MD results and continuum hydrodynamic predictions, even for the smallest systems undergoing 
single file transport of water. An optimum of permeability is found for an opening angle around 5°, in 
agreement with continuum predictions, yielding a permeability five times larger than for a straight nanotube. 
Moreover, we find that the permeability of hourglass shape nanopores is even larger than single nanopores 
pierced in a molecular thin graphene sheet. This suggests that designing the geometry of nanopores may help 
considerably increasing the macroscopic permeability of membranes. 


I. INTRODUCTION 

The study of flow through nanoscale channels has gen¬ 
erated a large activity over the last decades 1 ®. On the 
experimental side, transport across biological channels, 
as well as through artificial nanopores such as solid- 
state nanopores, nano-slits, or carbon/boro n-n itride nan¬ 
otubes, have been thoroughly investigated®. This has 
unveiled a number of new behaviors, such as the fast flows 
through carbon nanotube membraneiP®, non-linear elec- 
trokinetic transport^ or large osmotic transport in 
boron nitride tubefP. On the theoretical side, the ques¬ 
tion of fluidic transport at the nanoscale has generated 
a considerable activity, involving in particular molecu¬ 
lar dynamics (MD) simulations of flows in nanochan¬ 
nels. Since the pioneering work of Hummer^, nanochan¬ 
nels made of carbon nanostructures have been specif¬ 
ically explored- 11 . This has led to the recent high¬ 
lighting of the special properties of flows inside carbon 
nanotubesJ-Ldd^J and more recently through nanopores 
pierced in graphene 15 16 . 

Beyond the sole transport properties in nanoconfine¬ 
ment, the question of entrance effects is of particular 
importance for pores with nanoscale dimensions. This 
was discussed exhaustively in the phys iology literature 
since the work of Hille and Hall 17 *I2J, in the context 
of the entrance contribution to the electric resistance 
in ionic channels. However, hydrodynamic entrance ef¬ 
fects for the flow resistance into a pore were already dis¬ 
cussed a century ago by Sampson who calculated the 
flow resistance across a circular pore in an infinitely 
thin membrane, within the framework of continuum 
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hydrodynamief 19 . It was then generalized to a circular 
cylinder with finite lengtlP^ 21 . 

Hydrodynamic entrance effects involve a supplemen¬ 
tary viscous dissipation, which is expected to be domi¬ 
nant and limitating in systems of small length as com¬ 
pared to their radius, and/or systems with low friction 
(high slippage) at their walls. In this context, we showed 
in a recent contributiorP^l, using finite element (FE) cal¬ 
culations and theoretical estimates, that the hourglass 
geometry with a small cone angle reduces considerably 
the entrance dissipation and accordingly increases the 
permeability of the nanopore. Interestingly, the corre¬ 
sponding geometry matches quite nicely the shape of 
aquaporins, a biological water filter known to exhibit 
large flow rate^HSsl Now, biological nanopores, such 
as aquaporins, usually involve a single-file transport of 
water inside the core of the channel and continuum hy¬ 
drodynamics is accordingly e xpecte d to fail in the core of 
such sub-nanometric systems^ 2 This raises therefore 
the question of entrance effects at the interface between 
macroscopic reservoirs and single-file systems. How do 
entrance effects behave in such geometries? Does the 
hourglass shape retains its enhanced transport behavior, 
as suggested by continuum calculations? 

In the present paper we address these questions by 
exploring water transport across artificial hourglass-like 
channels using MD simulations. We show that - in a 
quite unexpected way - the results for transport in molec¬ 
ular pores compare quantitatively with continuum hydro¬ 
dynamics predictions. In particular, an optimized per¬ 
meability is still found for a small cone angle around 5 °. 
We then compare the hydrodynamic permeability across 
hourglass shaped pores to that of pierced graphene and 
show that the hourglass shape is significantly more effi¬ 
cient in terms of water transport. 












II. HYDRODYNAMIC RESISTANCE, ENTRANCE 
EFFECTS AND SAMPSON FORMULA 
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In the framework of weakly out-of-equilibrium sys¬ 
tems, the hydrodynamic resistance 7?hyd characterizes 
flow transport across a given structure, and is defined 
as the ratio between the pressure drop A P and the cor¬ 
responding flow rate Q between two reservoirs connected 
by the structure: 

-Rhyd = _ q ~ ■ (1) 

For a cylindrical channel with radius a and length L, the 
“inner” hydrodynamic resistance, relating the pressure 
drop along the channel to the flow rate, is given within 
continuum hydrodynamics by the well-known Poiseuille 
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where ry is the liquid viscosity, and a no-slip boundary 
condition (BC) has been assumed at the channel wall. 
In the case of liquid/solid slip at the pore wall, dissi¬ 
pation and accordingly the hydrodynamic resistance are 
reduced^: 
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where b is the so-called slip length 2 ! 

Beside this inner resistance, entering the pore gener¬ 
ates a supplementary hydrodynamic resistance, originat¬ 
ing from the bending of the streamlines towards the pore. 
This question was first discussed by Sampson in an ar¬ 
ticle in 1891^, where he calculated the velocity profile 
of a liquid flowing through an infinitely thin membrane 
pierced with a circular hole, using continuum hydrody¬ 
namics. He obtained the relationship between the flow 
rate Q and the pressure drop A p as 

Q = p, (4) 


where a is the hole radius and (7 = 3. Later, it was shown 
that Sampson’s formula could describe quantitatively the 
access pressure drop (i.e. the excess pressure drop due 
to the entrances) through a cylindrical pore, assuming 
the standard hydrodynamic no-slip BC of the liquid at 
the pore walls 2 ™ Therefore the access resistance to a 
cylindrical pore is given by a formula a la Sampson, as 
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with C a numerical constant, C ss 3. 

Gathering contributions, the total hydrodynamic re¬ 
sistance of a cylindrical pore with finite length is given 
as: 
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FIG. 1. Snapshots of the systems considered. Up: Armchair 
(6,6) tube. Down: Biconical system [armchair (5,5) tube, 
Tc/tteff = 20 and a c = 3.38 A; see text for definition of the pa¬ 
rameters]. Oxygen atoms are colored in red, hydrogen atoms 
in white, and carbon atoms in gray. 


In the no-slip case, the entrance effect becomes domi¬ 
nant when the aspect ratio L/a is below a threshold value 
[. L/a\ c = Cn/8 ~ 1. In the large slip limit b —► oo, R[ n 
and R ou t share the same scaling with the pore radius. 
Therefore, independently of the pore radius , the access 
resistance become the limiting factor for a critical pore 
length L c such that L c = ^R-b. For L<L C1 R to t ~ -Rout 
and the pressure drop will be concentrated at the inlet 
and outlet, with a constant pressure along the pore, as 
observed in previous AID work! 31 ! 

Now these different estimates rely on the validity of 
the continuum hydrodynamics, which is expected to fail 
whenever the diameter of the pore is in the nanome¬ 
ter rangfM2ZI. Therefore we now explore the validity of 
the continuum predictions for entrance effects using MD 
simulations of flow across pores with sizes varying from 
nanometers down to a few angstroms. 


III. MOLECULAR DYNAMICS SIMULATIONS AND 
FINITE-ELEMENT CALCULATIONS 

A. Molecular dynamics 


In order to quantify entrance effects, we use geometries 
depicted in Fig. [I] two water reservoirs are separated by 
a membrane consisting of two parallel graphene sheets, 
and pierced with a nanopore. The pore can be a sim¬ 
ple carbon nanotube (radius a c , defined by the center 
of carbon atoms, and length L = 10 a c ) or a biconical 
nanochannel, with a central nanotube and two conical 
entrances, made of graphene-likc sheets. The inner tube 
radius a c is varied from 25 A down to 3.5 A. We con¬ 
sidered both armchair and zigzag nanotubes. The po¬ 
sitions of the carbon atoms of the membrane are fixed 
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(simulations with flexible and fixed walls were shown to 
give si milar re sults for the statics and friction of confined 
liquid J 27 ^ — l ) . We use two pistons made of graphene 
sheets in order to maintain each reservoir at the desired 
pressure. Nanotubes are made by rolling a graphene 
sheet with respect to unit-cell. Periodic boundary con¬ 
ditions are imposed in all directions. The box size along 
the z direction is taken sufficiently large to ensure that 
the periodic images do not interact. In order to avoid 
hydrodynamic interactions between the tube and its pe¬ 
riodic images in the plane of the membrane, we choose 
a box with lateral dimensions (a :,y) equal to 10 times 
a c . We also make sure that reservoirs are bigger than 10 
times a c along z. There are 3.2 k water molecules and 
2.7 k carbon atoms for the smallest system ( a c = 3.5 A), 
and 800 k water molecules and 78 k carbon atoms for the 
biggest system (a c = 25 A). Finite element calculations 
indicate that, in that configuration, the error due to fi¬ 
nite size effects should be lower than 0.25 % (see following 
section). 

The Amber96 force fielcP^was used, with TIP3F® 3 wa¬ 
ter and water-carbon interaction modeled by a Lennard- 
Jones potential between oxygen and carbon atoms, with 
parameters eoc = 0.114 kcal/mol and croc = 3.28 A. 
There is no need to define a potential between car¬ 
bon atoms since they are fixed. The simulations were 
performed using LAMMPS^. Long-range Coulombic 
interactions were computed using the particle-particle 
particle-mesh (PPPM) method 3 ' 38 ( and water molecules 
were held rigid using the SHAKE algorithnP3. The equa¬ 
tions of motion were solved using the velocity Verlet al¬ 
gorithm with a timestep of 2 fs. 

Water molecules were kept at a temperature of 
300 K using a dissipative particle dynamics (DPD) 
thermostalpH This amounts to adding pairwise interac¬ 
tions between atoms, with a dissipative force depending 
on the relative velocity between each pair and a random 
force with a Gaussian statistics. This method has the ad¬ 
vantage of preserving hydrodynamics, even for complex 
3-dimensional flows as the ones considered here. The am¬ 
plitude of the dissipative term was carefully tuned to en¬ 
sure that the liquid viscosity is negligibly affected by the 
thermostat (this thermostating method was extensively 
tested and compared to other approaches in a previous 



Water molecules are initially disposed on a simple cu¬ 
bic lattice with equilibrium density. Pressure differences 
are imposed using the reservoir pistons, and a steady- 
state flow quickly appears, within less than 200 ps (see 
Fig. §. Note that this timescale matches the expected 
time for momentum diffusion in reservoirs, in agreement 
with the dominant role of entrance resistance. We then 
measure the steady-state flow rate by counting the num¬ 
ber of water molecules crossing the tube. This is shown in 
Fig. [2]for several tubes under a given pressure drop. The 
flux Q is deduced from the linear fit of the time depen¬ 
dent variation of the number of crossing water molecules 
ANu(t): Q = M/(pJ\fA_)dANii/dt , with M and p the 



t / ns 


FIG. 2. Number of crossing water molecules toward the right 
reservoir ANr as a function of time, for different tube radii 
o c (o c = 3.51 A in yellow, 4.31 A in cyan, and 4.69 A in 
purple. Full lines represent the MD results, and dashed lines 
represent the linear fit used to compute the flow rate. 

molar mass and density of water, and A/a the Avogadro 
constant. Finally, the hydrodynamic resistance is com¬ 
puted as R = AP/Q. For each geometry we ran a num¬ 
ber (up to 10) of independent simulations from different 
initial conditions, in order to estimate and reduce sta¬ 
tistical uncertainties. The production times ranged from 
~ 1 ns for the largest pores, to ~ 6 ns for the smallest 
ones. Although the results presented in this article were 
obtained for a pressure drop of 1000 atm — 1 atm, we em¬ 
phasize that linearity between flux and pressure drop has 
been checked systematically in all our simulations. 

Slippage effects - Independently from the above cal¬ 
culation, we also calculated the slip length of water on a 
graphene sheet. To this end, we performed simulations 
of Couette and Poiseuille flows of water confined be tween 
two graphene planes, in line with previous work 4142 . 
We measured a quite large slip length b = 123 ± 21 A. 
This means that tubes with nanometer range diameters 
can be considered as perfectly slipping. We also mea¬ 
sured the water viscosity, and found it to be equal to 
0.31 ± 0.02mPa.s, in good agreement with the expected 
value of 0.31mPa.s for the TIP3P model at 300 K®. 


B. Finite element calculations 

In order to compare MD results with continuum pre¬ 
dictions, we perform FE calculations using the commer¬ 
cial software COMSOL. We solve the Stokes equation 
r)Av = Vp in a 2D-axisymmetric geometry. For the 
cylindrical pore, the system is composed by a tube of 
radius a, and length L = 10a. Similar calculations are 
performed with the hourglass shape'll, see Fig. [3j The 
pore is connected to two reservoirs of characteristic size 
L r . Far from the pore, we impose a difference of pres¬ 
sure Ap between the two reservoirs and then measure 
the water flow across the tube, which gives us the hydro- 
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FIG. 3. Cylindrical and hourglass geometries used for the 
FE calculations. Two reservoirs of size L r are separated by 
a membrane containing a channel of radius o e fr and length 
L. The black dashed line corresponds to an axisymmetric 
condition. The full line represent the liquid/solid interface, 
we use perfect slip BC along it. We impose the pressure on 
the dotted lines. A c is the opening radius and a the cones 
angle. 



FIG. 4. Resistance of a single nanotube R = AP/Q, nor¬ 
malized by the bulk liquid viscosity r/, as a function of the 
radius o c of the tube. Circles represent MD results with arm¬ 
chair tubes, squares represent MD results with zigzag tubes. 
Lines are hydrodynamic predictions using FE calculations, for 
a c — a e ff = 3.2 A (dashed dotted line) and a c — a e g = 1.7 A 
(dashed line). 


dynamical resistance R = Ap/Q. Along membrane and 
tube walls we use perfect slip boundary conditions. Note 
that a small but finite friction may be expected at the 
water/carbon surface. However using typical slip lengths 
for the considered systems, as reported in the previous 
section, we checked that the corresponding effect on the 
hydrodynamic resistance is negligible as compared to en¬ 
trance effects: typically, the finite slip length affects the 
total resistance by ~ 0.2% for a c = 3 A and less than 
10% for a c = 22 A, and thus does not bring significant 
changes. 

Note that we use two different values of the radius a: 
a c , which is the distance between the centers of carbon 
atoms, and a e ft, which is the hydrodynamic radius of the 
pore, see Fig. [3] Typically one expects the hydrodynamic 
boundary condition (BC) to apply within the fluid at 
a distance of ~ oocl be. one molecular size from the 
walpS. Accordingly we expect for the radius a e g ~ a c — 
croc (this will indeed be in agreement with our results). 
However in the following we will merely use a e g as a 
free parameter, in order to identify the best agreement 
between the MD and FE results. For consistency, we 
introduce chamfers at the wall corners, with a radius of 
curvature a c — a e g. 

We also used FE calculations to estimate the minimal 
simulation box size needed to keep finite size effects to an 
acceptable level. Indeed, far from the entrances, the flow 
is identical to the one from a point source on a plane wall, 
implying that the liquid velocity decreases algebraically 
and size effects due to the finite size of the reservoirs in 
the MD simulations could therefore affect the flow rate 
versus pressure drop measurements. We checked that a 
characteristic size of reservoir L r ~ 10a o ff (a e g the pore 
radius) ensures an error of ca. 0.25 %. 


IV. ENTRANCE EFFECTS AT THE MOLECULAR 
SCALE: THE CYLINDRICAL PORE GEOMETRY 


Let us first explore the flow resistance across cylinders 
of finite length. Following the procedure described in 
section III A the hydrodynamic resistance R = A P/Q is 
extracted for the various tube diameters using MD sim¬ 
ulations. We plot in Fig. [4] the resistance normalized by 
the bulk viscosity, R/r /, as a function of the tube radius 
a c . These data are compared with continuum predic¬ 
tions obtained using FE calculations in the same geome¬ 
try. We used two different values for the hydrodynamic 
pore radius: one where the perfect slip hydrodynamic 
BC is located right at the surface of carbon atoms, i.e. 
a c — a e ff = <r c /2 ps 1.7 A; and the other where the perfect 
slip hydrodynamic BC is located at the first layer of wa¬ 
ter along the wall a c — a e g ps croc ~ 3.2 A, in line with 
results from Ref ! 41 * 44 ! As shown by Figure |4j we obtain 
a very good agreement between the MD results and the 
continuum hydrodynamic predictions obtained by FE. In 
line with previously quotted expectations, we found that 
a reasonable choice for the value of the hydrodynamic 
radius is a e fi — a c — 2.5 A. Note that nanotube’s chirality 
has no significant influence on entrance effects. 


It is particularly interesting to observe that the hy¬ 
drodynamic prediction is valid for nanotubes with effec¬ 
tive radii well below one nanometer, even when single file 
transport occurs. In such case, we expect a full break¬ 
down of hydrodynamics. However, entrance effects origi¬ 
nate mainly from the bending of the streamlines outside 
the tube, which occurs in the bulk on length scales much 
larger than the tube radius. This certainly explains the 
robustness of hydrodynamics to predict entrance effects. 




























V. ENTRANCE EFFECTS AT THE MOLECULAR 
SCALE: HOURGLASS NANO-PORES 
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A. Molecular flows through the hourglass geometry 


We now turn to the hourglass geometry. In our previ¬ 


ous work^, we showed on the basis of continuum hydro- 
dynamic calculations that the hourglass shape, inspired 
by aquaporins, lead to a large increase of the overall chan¬ 
nel permeability as compared to the cylindrical geometry. 
Here we explore how this conclusion can be extended to 
the case where the liquid in the middle part of the chan¬ 
nel is subjected to single file transport, as occurs e.g. in 
aquaporin channels. 

As discussed above, we build up a biconical nanochan¬ 
nel using graphene-like walls, see Fig. [TJ The cones ex¬ 
hibit an angle a, which is varied between 0 and 20 °. Con¬ 
ical entrance are made alike by rolling graphene sheet, 
thus leaving now a structural defect line. Note that the 
overall pore permeability is not affected by this defect 
line as tested using alternative methods to generate the 
cones. We choose a configuration in which there is single¬ 
file flow inside the central part, as found in aquaporins. 
Looking for possible departures from continuum hydro¬ 
dynamics, we consider the most extreme case, with the 
smallest radius that allows water molecules to fill the 
tube (a c = 3.38 A). Inspired by aquaporins, we choose 
cones’ lengths L c « 20 x a e ff = 17 A. We use reservoirs 
with a size at least ten times the opening radius, A c (with 
A c = a e ff + L c tana), which is sufficient to reduce finite 
size effects down to a negligible level, as estimated by FE 
calculations (see section III B). 

The MD results are shown in Fig. [5] We find that 
the hydrodynamic resistance is again minimized in the 
hourglass geometry for a cone angle a ~ 5°, even if the 
transport in the center part of channel is single-file. The 
minimal resistance is ca. 5 times smaller than the one of 
a tube with a straight entrance (a = 0). Furthermore, 
as demonstrated in Fig. [5j the continuum hydrodynamic 
predictions exhibit a good agreement with the results of 
the MD simulations for the hydrodynamic resistance. In 
particular, we highlight the critical importance of the po¬ 
sition of the hydrodynamic BC by showing two different 
cases of hydrodynamic radius for the inner part of the 
pore. 

We finally compare the results with our theoretical pre¬ 
diction for the hourglass resistance obtained in RefP^l, 
which writes 
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We used the values L c /a e g = 20, and the parameter 
C = 1.15 obtained from previous FE calculations in the 
case of an aspect ratio a e ff/(a c — a e ff) = 0.33 (see Ap¬ 
pendix) . The effective radius of the biconical channel was 
fixed to a e ff = 1 A= a c — 2.4 A in order to get best agree¬ 
ment. With this value, Eq. (J7| reproduces quite well the 



FIG. 5. MD results for the hydrodynamic resistance R = 
A P/Q of an hourglass channel (normalized by the bulk vis¬ 
cosity g) as a function of the angle of the biconical vestibules. 
Circles represent MD results. The dotted line represents 
the prediction of Eq. (JT]), with a e g = lA. The lines show 
the FE calculations with two different hydrodynamic radii: 
a e ff = 1.08 A (dashed line), a e g = 0.68 A (dashed-dotted line). 


MD results (see Fig. [5]). Note that MD results are best 
fitted with FE calculations for the same value of a e g (not 
shown for clarity). Therefore, the “hydrodynamic size” 
of the wall atoms obtained with the biconical geometry, 
a c — a e ff = 2.4 A, is quite close to the one obtained with 
a cylindrical pore in previous section (a c — a e g = 2.5 A). 
Additionally, one may notice that a number of MD re¬ 
sults are slightly outlying the main tendency, beyond the 
error bar. We conjecture that these may be attributed 
to discrete effects which may occur preferably for specific 
values of the cone angles. 


B. Comparison with pierced graphene 

Finally, we compare the efficiency in terms of trans¬ 
port, of hourglass shape nanopores versus pierced 
graphene. Pierced graphene was discussed lately to 
represent a very efficient geometry for desalination 
purpos (FE6l45l , 

in particular due to its large permeabil¬ 
ity combined with a excellent rejection ability. The large 
permeability of pierced graphene is due to the molecular 
thickness of the graphene sheet: in this case, the trans¬ 
port is fully controlled by the entrance effect, and belongs 
to the Sampson’s class of problem. It is therefore inter¬ 
esting to compare the performance of this geometry to 
the hourglass one which was precisely found to reduce 
entrance effects. 

To this end we have performed MD simulations of 
transport across nanopores drilled in a graphene sheet 
with various pore sizes, and measured the correspond¬ 
ing hydrodynamic resistance. This is compared to the 
hydrodynamic resistance of an hourglass shape nanopore 
with an inner channel having the same radius. We use 
hourglass systems with three different inner pore radii, 
respectively a c = 3.38, 4.06 and 4.73 A. We keep the 
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FIG. 6. Comparison of hydrodynamic resisttance of hour¬ 
glass, biconical systems (purple triangles) and graphene 
sheets pierced with a single hole (orange circles). Dashed lines 
are continuum hydrodynamics predictions (see text). The 
cyan dashed line is the Sampson formula for a hole in a in¬ 
finitely thin membrane. 


cone lengths ratio L c /a c equal to 20 and the angle a 
equal to 5°. The radius of the hole a c in the graphene 
sheet is equal to 3.38, 4.06 and 4.73 A. Results are shown 
in Fig. [6j As can be seen in this figure, the hourglass 
nanopore has a hydrodynamic resistance which is sys¬ 
tematically smaller than the graphene for the same inner 
pore diameter. As a comparison, we plot the result of 
the classic Sampson formula, Eq. Q, with a hydrody¬ 
namic radius equal to a e g = a c — 2.5 A (see above) and 
a coefficient C = 3. As seen in Fig. [6j this overesti¬ 
mates the MD results for the hydrodynamic resistance of 
porous graphene. In an attempt to improve this result, 
we take into account at the continuum level the finite 
thickness of the graphene membrane, as compared to the 
pore radius (see appendix for details). This leads to a 
reduction of the coefficient C in the Sampson formula. 
The value of C depends on the ratio a e ff/(a c — a e ff), and 
in the present conditions it varies between 1.15 and 1.6 
for the considered radii. We choose C ~ 1.3 in Eq.Q as 
a compromise, showing a good agreement with the AID 
results. For the hourglass geometry, we compare the AID 
results with the prediction of Eq. ([7]), using the same 
coefficient C = 1.3. This shows again a good agreement 
with the AID results. 

Altogether these results show that tuning the geometry 
of nanopores allow to strongly optimize water transport 
through membranes. The hourglass shape outperforms 
both nanotubes and pierced, molecular thick, graphene. 


VI. CONCLUSION 

In this paper, we investigated the fate of hydrodynamic 
entrance effects in the limit of molecular scale nanopores, 
down to the single-file regime. Strikingly, we found 
an extremely good agreement between AID simulations 


and continuum hydrodynamics predictions of entrance 
resistance obtained from FE calculations. For straight 
nanopores, comparison between AID and FE approaches 
incorporate the position of the solid-liquid interface as an 
important parameter, which therefore could be located 
using such an approach. In addition, the agreement be¬ 
tween molecular scale transport and continuum hydrody¬ 
namics extends to more complex geometry of nanopores, 
inspired by the hourglass shape of aquaporin biological 
nanochannels. 

This hourglass shape proves to substantially improve 
the water transportation efficiency down to single-file 
regime, with an optimum for shallow opening angles 
~ 5°, all in line with continuum hydrodynamics predic¬ 
tions. Compared with the promising system made of a 
circular nanopore drilled in graphene sheets, or straight 
carbon nanotubes, the hourglass shape pore appears far 
more efficient, illustrating the dominant role of entrance 
dissipation in all these systems. Overall, this stresses the 
necessity for finely tune the geometry of nanopores for 
strongly optimizing water transport across membranes, 
a task for which simple continuum approaches can be 
astonishingly reliable. 
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APPENDIX: SAMPSON FORMULA AND MEMBRANE 
SIZE EFFECT 

Carbon nanotubes are made of atoms with a finite ra¬ 
dius, and the graphene membrane has therefore a finite 
molecular width which could compare to the diameter of 
the pore, for nanometric pores. At the level of continuum 
description, this finite aspect ratio is expected to modify 
the prefactor C in the Sampson equation, Eq. ©■ We 
estimate here the modified coefficient C using FE cal¬ 
culations. The effective radius of carbon atoms (defining 
the hydrodynamic radius of the pore) is equal to a c — a e g, 
see figure [3] in the main text. We explored a large range 
of radius a e g for fixed a c — a e g. Figure [7] shows the effect 
of the aspect ratio on the value of the C coefficient, for 2 
different membrane thickness. As we can see, for a ratio 
fleft/(flc~ a eff) ~ 1, which correspond to a (5,5) armchair 
tube, C is around 1.5 instead of 3.75, corresponding to 
a 60 % improvement of the pore permeability. Note that 
this upper value of C = 3.75 corresponds to the thick 
membrane limit achieved for our hourglass pore and that 
Sampson’s prediction C = 3 is recovered for vanishingly 
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FIG. 7. FE testing of the Sampson formula with a mem¬ 
brane of finite thickness: Measured Sampson coefficient as 
a function of the ratio effective radius a e ff over chamfer ra¬ 
dius a c — a e fr- (long dashed line) nanotube limit L/a e s = 5; 
(dashed line) atomically thin membrane nanopore limit L = 
0 . 


thin membranes as can be seen in figure [7] for L = 0. 
Surprisingly, the frictionless inner pore section increases 
the global, entrance-associated, dissipation as discussed 
in our recent work 22 . 
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